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Abstract 

Starting from the T-Q equation of an open integrable spin-^ XXZ quantum spin 
chain with nondiagonal boundary terms, we derive a nonlinear integral equation (NLIE) 
of the sine-Gordon model on a finite interval. We compute the boundary energy and 
the Casimir energy for the sine-Gordon model with both left and right boundaries. 
A relation between the boundary parameters of the continuum model and the lattice 
model is given. We also present numerical results for the effective central charge of an 
open spin-i XXZ quantum spin chain which find agreement with our analytical result 
for the central charge of the sine-Gordon model in the ultraviolet (UV) limit. 
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1 Introduction 



Due to applications in statistical mechanics and condensed matter physics, spin-^ XXZ 
quantum spin chain and sine-Gordon models with two boundaries have been subjected to 
intensive studies over the years [l]-[20]. Works on these topics have covered mainly diagonal 
and up to certain extent, nondiagonal boundary interactions. Proposal of Bethe ansatz 
solutions for the XXZ quantum spin chain with nondiagonal boundary terms [T71 dB] have 
made it possible to study these models further [191 ISO]- One such study is the investigation 
of the finite-size effects in sine-Gordon models using the nonlinear integral equation (NLIE) 
approach [211 122]. Also refer to [9l [23] for further work on the subject. 

Motivated by previous works on the subject, we use a recently proposed solution of an 
open spin- 1 XXZ quantum spin chain [25] to compute finite-size effects to the lowest energy 
state of the sine-Gordon model with two boundaries. In contrast to previous works, our 
work involves a boundary parameter free of any constraint (describing nondiagonal boundary 
terms). However, one limitation of this solution is that the bulk anisotropy parameter for 
the XXZ spin chain assumes special values, namely ifi, where = ^, z/ = 3, 5, 7, . . .. In 
particular, we derive a NLIE for the lowest energy state of the sine-Gordon model on a finite 
interval and compute the boundary and Casimir energies. Our motivation is two fold: While 
the lowest energy state for other studied (critical) spin-^ XXZ quantum spin chain models is 
described by a sea of real Bethe roots, the lowest energy state for the spin-| XXZ chain model 
considered here is described by a sea of "two-strings", i.e., complex conjugate pairs of Bethe 
roots, which is rather a characteristic of spin-1 XXZ chain. We thus feel it is worthwhile 
and interesting to compute finite size corrections for such a spin-^ model. Moreover, such a 
computation will serve as a useful guide when one considers the corresponding open spin-1 
XXZ chain which has been associated with supersymmetric sine-Gordon models [25j-[31]. 
As pointed out in [21], due to the sea of "two-strings", familiar method of deriving the NLIE 
based on Bethe ansatz equations and the counting function [22] does not seem to work. 
Fortunately, an NLIE can still be derived from the model's T — Q equation. Refer to [32] 
for more details. 

The outline of this article is as follows. In section 2, we briefly review the sine-Gordon 
model and the Hamiltonian of the open spin-i XXZ quantum spin chain with two boundaries 
[T9l [20]. In section 3, as a warm up exercise, we rederive the NLIE given in [19] for the open 
spin-i XXZ /sine- Gordon model where the boundary parameters obey a certain constraint. 
However, in contrast to the approach [22] used in [191 120] , we employ a method that utilizes 
the T — Q equation [32] of the open spin-| XXZ quantum spin chain model. In section 4, 
we give main results of the paper. We derive an NLIE for the sine-Gordon model again 
using the T — Q equation for an open spin-i XXZ quantum spin chain with nondiagonal 
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boundary terms, but now one of the boundary parameters free of any constraint [23]. We 
also derive the boundary energy and Casimir energy for the lowest energy state for this case. 
Finally, we show that the analytical result derived for the central charge of the sine-Gordon 
model in the ultraviolet (UV) limit agrees with the numerical results for the central charge 
of an open spin-i XXZ spin chain that we obtain by solving the XXZ spin chain model 
numerically for few finite number of sites A^, and extrapolating the results to — )■ oo limit 
using an algorithm due to Vanden Broeck and Schwartz [33]- [35]. Also see [2], where such 
an extrapolation technique was used to study certain properties of statistical mechanical 
systems. This is followed by a brief discussion of our results and some open problems in 
section 5. 



2 The sine-Gordon model and open spin- 1/2 XXZ spin 
chain with two boundaries 



In this section, we briefly review the sine-Gordon model (reproduced from [T9|) and the 
Hamiltonian of the open spin-i XXZ quantum spin chain with two boundaries. The sine- 
Gordon model on the finite "spatial" interval x G [x_ , a;+] is described by the (Euclidean) 
action 



dy / dx A{(p , df,(p) + / 

-oo J X— J —c 



where the bulk action is given by 



dy 



dip, 
dy' 



dip, 
dy' 



A{ip , df,ip) = -{d^ipf + Hbulk COs(/3(y9) 



and the boundary action is given by 

dip. 



B±iip^ 



dy 



l.^oos{-iip-^,))±^- 



(2.1) 



(2.2) 



(2.3) 



As noted in [191 [20], the action is similar to the one considered by Ghoshal and Zamolodchikov 
[6], except for two boundaries instead of one. Moreover, the presence of an additional term 
depending on the "time" derivative of the field in the boundary action (12. 3p is also well 
noted. As mentioned in [12], while such a term can be eliminated in the one-boundary 
case by including a term proportional to d^dyip in (12. 2p . such a step would simply result in 
the elimination of only one of the two 7-1- parameters (say, 7+), and in a shift of the other 
(7_ I— 7_ — 7+) in the two-boundary case. The continuum bulk coupling constant (3 is 
related to the lattice bulk coupling constant /i by = 8(7r — /x) = 87r(z/ — l)/z/, taking 
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In subsequent sections that follow, we shall consider the energy of the ground state of this 
model as a function of the interval length L = x+ — a;_, for large L. The leading contribution 
which is of order L and which does not depend on the boundary interactions is well known 
[36] . The boundary correction of order 1 is also known pT] [T2j . In this paper, we shall 
mainly concentrate on the compution of the Casimir correction of order and derive a 
nonlinear integral equation [2T]-[23] for the lowest energy state. The length L and the soliton 
mass m (whose relation to ^ibuik is known [36J) are given by 

L = NA, m = ^e'^, (2.4) 

respectively. In (12. 4p . A is the lattice spacing, which in the continuum limit, taken to be 
A —7- together with A — t- oo and N ^ oo for the inhomogeneity parameter and number 
of lattice sites respectively. We shall return to this in following sections where the NLIE 
are derived. Further, as given in p2], the boundary parameters in the continuum action 
(/i± , v^o ) is related to the boundary parameters of the lattice model {a± , (3±) (that appears 
in the Hamiltonian of the open spin-^ XXZ quantum spin chain (see (12. 7p below)) by 

sinh(a± + /3±) = i^te-i^^o ^ 



sinh(a±-/3±) = ^te+^^^o , (2.5) 



where 



(J'C 



2/i, 



bulk 



(2.6) 



sin ^ 

The Hamiltonian of the open spin-i XXZ quantum spin chain is given by O [6] 

^ N-1 

= 2 { («+l + ^n^n+l + cosh T] 

n=l 

+ sinh 7] ^ coth a_ tanh /S-crj" + csch a_ sech /3_ ( cosh 9-ai + i sinh O-cf) 

— coth tanh + csch a+ sech /3+ ( cosh 9+a% + i sinh O^a^) | , (2.7) 

where , are the usual Pauli matrices, 77 is the bulk anisotropy parameter, a± , /3± , 6± 

are the boundary parameters, and N is the number of spins. 

We remark that (12.51) was derived for the case where the lattice parameters {a± , P±,6±) 
obey the following constraint 

a- + ^- + a+ + /3+ = ±{e--e+)+r]k, (2.8) 
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where k is an even integer if N is odd, and is an odd integer if is even. A convenient 
redefinition of bulk and boundary parameters can also be adopted [19] : 

7] = ifi , a± = ifia± , I3± = fib± , 9± = i^c± , (2.9) 

where /i , a± , fe± , c-t are all real, with < < vr. With the above redefinitions, the constraint 
relation (12. 8p yields the following pair of real contraints: 

a_ + a+ = ±|c_ — c_|_| + , 

h- + h+ = 0. (2.10) 

In this paper, we shall consider only even case. We also remark (see [19] for details) that 
the lattice parameters 6± are related to the continuum parameters j^. which appear in fl2.3p . 
We shall see that for the case studied in section 4, (12. 5p still holds true. 



3 NLIE of the spin-| XXZ / sine-Gordon with constraint 
nondiagonal boundary terms 

In this section, we shall rederive the NLIE for the spin-i XXZ /sine- Gordon model (along 
with the order 1 and order 1/L correction to the energy) with constraint (12. Sp among the 
lattice boundary parameters. In contrast to the familiar approach [22] used in [19], we utilize 
instead the model's T — Q equation (describing the transfer-matrix eigenvalues T{u)) and 
its analyticity properties. We shall follow closely steps utilized in [31] and employ similar 
notations. 



3.1 T -Q equation and NLIE 

The T — Q equation of the inhomogeneous open spin-| XXZ chain with general boundary 
conditions (but with the boundary parameters (a± , f3±,0±) obeying constraint (12.80 ) is given 
by [18] 

T{u) = sinh(2w + 2/x)5W(«)0(w+^)%f^ 

2 Q{u) 

+ smh{2u - z/i) B^-\u) <p{u - ^) ^(;; + ^/^) , (3.1) 

2 Q[u) 
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where 



(f){u) = sinh'"^ (m - A) sinh'' (m + A) , 

'[u) = smh[u± — - — ) smn(u ± — - — ) cosn(u =F — - — j cosn(u =F — 7: — J, 

^ Zi 

M 

= W sinh('U — ffc) sinh(-u + f fc) . (3.2) 
fc=i 

where the bulk anisotropy parameter is 77 = i/x, and we have redefined the boundary param- 
eters as A± = 2a± — 1 , B± = 2ih± + 1. A is the inhomogeneity parameter which provides a 
mass scale (see (12 ■4p ). Vk represents the Bethe roots which are also zeros of Q{u). We note 
here that the Q{u) above differs from that given in [18] by a mere shift of ^ but otherwise 
equivalent. 

Next, following [31j, we define the auxiliary functions a{u) and a{u) by 

sinh(2n + z^) B^ju) 4>{u + f ) Q{u - z/i) 1 

a[u) = , a[u) = a[—u) = , , . (3.3) 

sinh(2M-z/i)5(-)(M)0(M-f)g(M + z/i) a{u) 

The transfer-matrix eigenvalues then simply become 

T{u) = sinh(2u-z/x)5(-)(w)0(u-^)^^^^(l + a(u)) 

= sinh(2u + tfi) B^^\u) (P{u + !^) ^(^-^^) (1 + a{u)) , (3.4) 

Note that T(m) does not have zeros near the real axis except for one simple zero on real axis 
at M = 0. The Bethe Ansatz equations can be written as [18] 



a{vk) = -l. k = l,...,M. (3.5) 

We consider the lowest energy state with ^ real Bethe roots, namely M = We shall 
only consider "massless regime" (with purely imaginary bulk anisotropy parameter), rj = ifi, 
with < /i < TT. The regions in parameter space A± which yield real Bethe roots for the 
lowest energy state can be divided in the following way. (See [19] for further details on some 
numerical results about these parameter regions.) 



/ 
// 
/// 
IV 



0<A±<^ 

0<A, <^&-^<A_<-l 
, " 3.6 

2" < A± < -1 ^ ^ 



2^ < < -1 & < A_ < ^ 
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In addition, due to fl2.10p . the parameters B± satisfy 

B,+B_ = 2 



(3.7) 



Also as performed in [31], one can remove the root of T{u) at the origin (as pointed out 
earher, T{u) does not have zeros near the real axis except for a simple zero at the origin) by 
defining 

f{u)=^^, (3.8) 

where fi{u) is any function whose only real root is a simple zero at the origin, that is 
fi{0) = , /i'(0) 7^ 0, where the prime denotes differentiation with respect to u. Hence, the 
new T — Q equation becomes 



Q[u) Q[u) 



where 



t±{u) 



sinh(2M ± ill) ~ 



(3.9) 



(3.10) 



Utilizing the analyticity of InT(M) near the real axis, we have the following from Cauchy's 
theorem, 

= ^ rfw [lnt(u)]V'=", (3.11) 
where the contour C is chosen as in Figure 1, e being small and positive. 



Co 



Figure 1: Integration contour 
As a result, after using (13. 9p . (13. lip can thus be written as 

Q{u — ifi) 



= du [lnt+(M)]"e^^"+ / du {In 
Jci Jci 

+ f du [lnt_(u)]"e^'="+ f du {\n 
J C2 J C2 




e*'^"+ / du [ln(l + a(M))] e 

Ci 



" Aku 



e*'''"+ / du [ln(l + a(M))]"e 
C2 



" iku 



Following [31], we define Fourier transforms along C2 and Ci as 

Cr'{k)= [ du [ln/(w)]V'=" 



Lf"{k)= / du[\nf{u)]"e 
'C2 



II Aku 



(3.12) 



(3.13) 



Ci 
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respectively, and after exploiting the periodicity 

Q{u) = Q{u — in), u G Ci, and Q{u + ifi) = Q{u + ifi — iir), u E C2 
we have the following 

Q{u - ifi) 



(3.14) 



du < In 



Ci 



du < In 



C2 




Aku 



Aku 



-LQ"{k) (e-^*^ - e-^^) 
LQ"{k) [e^^'-^^^ - l] , 



Further, defining 



'Ci JC2 



" Aku 



we obtain from fl3.12p 

C{k) + LQ"{k) [e^^'-^^^ - 1 - e'^"^ + e-"'^] + LA"{k) + £l''(A;) = 
which leads to the following 



LQ"{k) 



e 2 



LA"(A;) + CA"{k) + C(A;) 



4cosh(f )sinh((7r-^)|) 
Note that the following definitions have been adopted, 

A{u) = 1 + a{u) , A{u) = 1 + a(M) 



Having found (IS.lSp . we proceed to derive the NLIE for the lattice sine-Gordon 
Taking the Fourier transform of a{u) (see (13. 3p ) along C2, one arrives at 



(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 
model. 



La"{k) 



+ 



1* \ 

















iku 



Q{u — ijj) 
Q{u + i^) 

sinh(2M + z/i) B^^'^{u) (j){u + 



sinh(2M - ifi) B(-){u) 0(m - f ) 



Aku 



(3.20) 



Using (13.141) . we obtain 



+ I du <.\n 



sinh(2u + ifi) B^+\u) (f){u + f ) 
sinh(2M - ifi) 5(-)(n) 0(u - f ) 



Aku 



(3.21) 



^This is to make the imaginary part of the argument negative. 
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Inserting fl3.18p into fl3.2ip yields the required NLIE for the sine-Gordon model in Fourier 
space, 

L^'{k) = G{k) \LA"{k) + CA"{k)] + Crik) , 



where 



G{k) 



sinh ((tt - 2/i)|) 
2cosh(^f;)sinh((7r-/i)|) 



Crik) 



G{k)C{k)+ I du {In 



sinh(2M + i/i) B(+\u) 0(m + f ) 
sinh(2u - ifi) B(-){u) 0(m - f ) 



(3.22) 



(3.23) 



e*'=" . (3.24) 



The second term in (13.241) yields 



du I In 



C2 



sinh(2M + 2/i) B^+\u) (j){u + f ) 
sinh(2M - ifi) B(-){u) (t){u - f ) 



27r?/'(fc)< s_e 



-m|a_| 



iku 



s+e 



Ak I —i\k\ 



+ A^(e*'"' + e 

where s± = sgn(yl±), 'ipik) 
identities (see also [3T]). 



e 2 



l-e 



fc and 4'2ik) 



+ 27r^2(A;) 



(3.25) 



and we have used the following 



[ — [In sinh(M — ia)]" e 
Jc2 



" Aku 



-k{a—nn) 



(3.26) 



where n is an integer such that < 3?e(a — nir) < tt, and 



C2 



nil 

— [lnsinh(2u)]"e''^" = ^k) 

271 



(3.27) 



Further, from f l3.10p . fl3.16p and using the fact that B^{u) and (piu ± y) are analytic and 
nonzero near the real axis (hence changing the contour integral on Ci into — C2), one arrives 
at the following result for C{k), 



C{k) 



— / du {hi 

'C2 



sinh(2M + ifx) S(+)(m) 0(m + f ) 
sinh(2M - ijj) B^-){u) (t){u - f ) 



Aku 



j) du [In /u(-u)]" e 



" iku 



(3.28) 



Note that the first term in (I3.28P is simply the negative of ( I3.25p . The second term in ( I3.28P 
reduces to —27ik. (Refer to our earlier discussion on /i(u)). Using fl3.23p -( l3^25l) and f l3.28p . 
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we finally obtain the following for Cxik), 

jNcos{Ak) s+sinh((/x|A+| -7r)|) +s_sinh((/i|A_| -7r)|) 

Gt fc = — Z7rfc< ; i ; — 

I cosli(f) 2cosli(f )sinli(( 

[sinh + sinh cosh(^) sinh ((2/i - 7r)|) 



(3.29) 



2 cosh(f ) sinh {{fi - 7r)|) cosh(f ) sinh ((/i - 7r)f ) 
Converting ( 13. 22^ to coordinate space and integrating twice, we obtain 

/oo roo 
du' G{u — + ie) ln(l + a{u' — ie)) — / du' G{u — u' — ie) ln(l + a{u' + ie)) 
-oo J — oo 

/ sinh ^ \ 

- z2Artan~M — ^ j +^PM,^(M) + Z7r, (3.30) 

In O3.30p . G{u) is the Fourier transform of G{k) (see (13.231) ) 

1 

G{u) = — dk e-^'^" Gik) , (3.31) 



2vr „ 

and Pbdryiu) is given by 

Pbdry{u)= du'R{u') = - du'R{u'), (3.32) 

Jo ^ J ~u 

where R{u) refers to the Fourier transform of R{k) which is given below, 

s+ sinh ((/u|A+| - 7r)|) + s_ sinh ((Ai|yl_| - 7r)|) 



R{k) = -2tt 



2cosh(^)sinh((^-7r)|) 
[sinh + sinh cosh(^) sinh ((2/i - 7r)|) 



2 cosh(^) sinh ((/i - 7r)|) cosh(^) sinh ((/i - 7r)|) 



(3.33) 



The factor ivr (integration constant) in f l3.30p is obtained by considering the m — )■ oo limit 
of (13. 3p and (I3.30p . Further, proceeding as in [31], one obtains the correct factor. We have 
explicitly checked that this procedure yields the same integration constant for all possible 
combinations of the boundary parameters, namely all four regions given in (13. 6p . 



sinh 2^ 
M 



Next, taking the continuum limit (A— >oo,A^— 7'00,A— i-O), the term —i2N tan ^ 
becomes —i2mL sinh 6 after defining the renormalized rapidity 6 as 

e = ^. (3.34) 
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Thus, fl3.3Up becomes 

/OO POO 
de' G{e -e' + te) ln(l + a{e' - te)) - / de' G{e -e'~ ie) ln(l + 3(6' + ie)) 
OO J —OO 



— i2mL sinh 6 + i Pbdryi^) + ^^r , 
where following definitions have been used, 

e = -, a{e) = a{^), P,,,y(e) = Puryi—) , G(0) = ^G(^) 
/i vr vr 7r TT 

By making following identifications 

f{e) ^ H-m) , nZ(^) = -p....(^) 



where Pbd^y{0) is labelled as Pbdry{(^) in tlBj, and defining instead, 



(3.35) 



(3.36) 



(3.37) 



9 = nu, G{9) = G{- 



1T 



(3.38) 



one arrives at the result found in [19] (see equation (3.26) of the reference), hence confirming 
that the result fl3.35p matches that of Ahn and Nepomechie. 

3.2 Boundary and Casimir energies 

In this section, we compute the boundary correction (order 1) and the Casimir correction 
(order l/L). We begin by following the prescription of Reshetikhin and Saleur [37] (see also 
[3T]). according to which the energy for the inhomogeneous case (A ^ 0) is given by 



E 



A\du ^ ^ 



u=A+f 



-^lnT(M) 

du 



(3.39) 



where g is given by 



Using the fact that 



— InT(M) 
du 



d 



u=A±i 



— InTW(M) 
du 



(3.40) 



(3.41) 



u=A±i 



where 



ijj, Q{u =F i/i) 
Q{u) 



(3.42) 
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and 



[In /(«)]'= / -Lf'{k)e-^^- 
Zn 



u e C2 



(3.43) 



(which in fact follows from f l3.13p ) f l3.39p reduces to 
E = 



4-^^lnTW(« + ^)-lnT( '(u 
Adu 



u=A 



9_ f dk ikA 
A I 27r^ 



Using (13.421) . one finds the following 



'e'^LTM'{k)-e-'TLT(^'(k) 



(3.44) 



{—ik) 

- e-^Hi^Yik) - e-'^'U'ik) - e'^^^^P^ 

{—ik) 



(3.45) 



But first one need to determine the explicit form for LQ'(k) which needed to be substituted 
in ([335]). After using (l338|) . (13:281) and LQ'{k) = j^^LQ"{k), we obtain 



LQ'ik) 



4cosh(f )sinh((7r-/i)|) I ^ ' ^ ' ^ ' ^ ' 



+ 



L^'{k) - [- 



27r^2(fc) 



- 27rz 



This eventually leads to the following result for the energy, 

E = El + El + Ei/L 

where 
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Er 



27rA 



,^k 



dk e-'^'^ia.iihi^—) 2e-- cosh((/i - -)k)L(f)'{k) 



(3.46) 



(3.47) 



(3.48) 



^1 



27rA 



'ikA 



tanh( 



fik. 



2m 2'Ki cosh((/i — ^)k) 



sinh ^ 



+ 



+ 



sinh ^ 



)'(fc) + e'^L5(+)'(A;) 



(3.49) 
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Ei/i 



27rA 



dk e 



-ikA 



cosh 



LA'{k) + CA'{k) 



(3.50) 



As the subscripts suggest, equations f l3.48p . f l3.49p and f l3.50p refer to the bulk, boundary and 
Casimir energies respectively. Now, we shall evaluate each of these terms explicitly. Using 



L(p'{k) 



L(t)"{k) 



—ik 



and the identity ( I3.26p . we have the following. 



which upon substitution in f l3.48p . yields 
El 



^ r dk (1 + e-^^^'^: 



-ik) ' 

sinh^cosh((/i- f)A;) 
cosh ^ sinh ^ 



(3.51) 



(3.52) 



(3.53) 



Adopting the renormalization procedure [22] of discarding divergent terms (A-independent 
term), keeping only the finite terms, evaluating the remaining integral by closing the contour 
in the lower half plane and selecting only the contribution from the pole at /c = — ^ and 
using (12. 4p and (I3.40p in the process, we finally arrive at the result found in [22| 136], namely 



^ 1. 2 

El = -Lm cot — 
4 2/i 



(3.54) 



We now consider the boundary energy I KWf . Using results fl33TD (with (p B^^^) and 
( I3.26p . we have 

2iTT 



e-'TLB'^-y{k) + e^LB(+y{k) 



sinh 



Trfc 



cosh[^-s_(^ + ^A_)] 



+ cosh[— -s+(^ + ^A+)]+cosh(^i?_ 
+ cosh(y5+-— ^ 



kfi^ 
2 ' 



Substituting fl3.55p into fl3.49p . we obtain the boundary energy 



El 



m 
~2 



nu cos (^(vr - s_/i(l + A_))) 
1 + cot ^-^ ^ 



sm 



B- 



sm 



sm 



(3.55) 



(3.56) 



which is evaluated using the same contour as for the bulk energy. The symbol 
represents the terms with A_ — )• , B_ B^ , s_ — > s+. In terms of a± and b±, the above 
becomes 



El 



m 
~2 



TTU COS (^(tc — 2s_/ia_)) cosh (6„7r) 
1 + cot ^ ^ ^ ^ 



4 sin ^ sin ^ 



(3.57) 
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which is the expression found by Ahn and Nepomechie in [T9] . 

Finally, we consider the Casimir energy given by f l3.50p . As with the NLIE, passing to 
coordinate space and taking the continuum limit, we obtain 

= ft L * ^"^ ( cosh ;(A - '"<^ + + 

Further, ( ^^.^^ ^ ^e^^^^""-^') at A oo limit. Using ([MD, we have 

Ec = / c/M5>meM("+"Mn(l + a(M + ie)) (3.59) 

2/^ i-oo 

From fl3.34p and the first two definitions in f l3.36p and after some manipulation, Ec reduces 
to (Note also that we have used A(u) = A(— n) and '^mz = —'^mz in the process.), 

Ec = / c/^S^m sinh(^ + ie)ln(l + a(^ + ie)) (3.60) 

Invoking the identification f{9) = ln(— a(6')), the above becomes 

Ec = / d9 3m sinh(^ + ie) ln(l - e^^^+^'^) . (3.61) 

27r J-oo 

which is indeed the result derived in 



4 NLIE of the spin-| XXZ/sine-Gordon with nondiag- 
onal boundary terms 

In this section, we give the main results of the paper. We derive a NLIE for an open spin-| 
XXZ/sine-Gordon model with the lattice boundary parameters satisfying the following, 

a± = r] 13+ = p_ = p e^ = e+ = Q. (4.1) 

Note that the above parameters do not obey the constraint (12. 8p . Moreover, we take /3 to be 
arbitrary and real. The bulk anisotropy parameter rj, equals i/i, where /z = ^ , z/ = 3, 5, 7, . . .. 
As in section 3, we will also derive the order 1 and 1/L corrections to the energy. We 
shall again use the spin-| XXZ quantum spin chain model's T — Q equation [21]. We then 
consider the UV limit of the central charge for the sine-Gordon model. We also present some 
numerical results. 
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4.1 T -Q equation and NLIE 



The T—Q equation of the inhomogeneous open spin-| XXZ chain with nondiagonal boundary 
conditions, with the boundary parameters {a± , P± ,0±) specified by (14.11) is given by [2] 



T[u) = sinh(2M + z/x) sinh(M H ) sinh(-u ) B^^'{u) (p{u -\ ) 

+ sinh(2w - z/i) sinh(« - ^) sinh(n + ^) B^-\u) <P{u - ^)^^'' - ^/^ " 1)) 



Qiu) 



(4.2) 



where 



0(m) = sinh^'' (m - A) sinh'^ (u + A) , 
i3W(«) = sinh(i(«T^))sinh(i(«T^))cosh(i(«T^)) 



X 



COsh(^(MT^^)) 



M ^ -1 

Q{u) = JJsinh(-(M-^fc- — ))sinh(-(M + Wfe + — )) 



k=l 



In (iJ), 



u = 3,5,7,..., M = Ar + z/-l. (4.3) 



A± = l + 2ib± - z/ , 5± = 1 - 2i6± + u . (4.4) 



where b± is related to /3± (that appear in spin chain Hamiltonian (12.71) ) as in (12. 9p . Note 
that /9_ = /9+ = /9 imphes 6+ = 6_ = b. This in turn yields 

A+ = A^ = l + 2ib-u, B+ = = l- 2ib + p (4.5) 

{"i^fc + y) represent the Bethe roots (and the zeros of Q{u)). The N + v — 1 "shifted roots" 
{vi , . . . , vm+u-i} for the lowest energy state have the following structure, 

7/), + HE ■ h — ^ 2 — 

where {wkiW^j^^ I'^f''} all real and positive, {wfe} are the sea-roots while {w^"'^'' 
represent the extra-roots, which are not part of the sea-roots. As pointed out earlier, an 
interesting feature of this model is that the lowest energy state is described by y "strings" 
of length 2, in addition to pairs of "strings" of length 1. By choosing /9+ = we find 
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that wl = wl . Numerical studies indicate that foi N ^ oo, wl becomes large as well. It 
is also worth mentioning that from numerical results for the Bethe roots, Wj^^ > Wk for all 
possible values of j and k, recalling the fact that j = 1 , 2 , . . . , and k = 1,2,... , y- 

Next, as in section 3, we define the corresponding auxiliary functions a{u) and a{u) by 

sinh(2M + ifi) sinh(M + f ) sinh(M - ^) B^+\u) 0(m + f ) Q{u + ifi{u - 1)) 
QjIii] = — — 

sinh(2u - ifx) sinh(u - f ) sinh(u + f^) B(-){u) 0(m - f ) Q{u - - 1)) 

d{u) = a{-u) = ^. (4.7) 
a[u) 

Using (14. 7p . the T — Q equation (14.21) becomes 

T{u) = sinh(2u — ifi) sinh(u — —) sinh(u + ^--^) B^~\u) (f){u — — ) 

2 2 2 

Q[u) 

= sinh(2'u + ijj,) sinh('u + — ) sinh(-u —) B'''^\u) (f){u + — ) 

Q(u + iuiu — 1)) . , , 

X ^1 + aw , 4.8 

Q{u) 

As before, T{u) does not have zeros near the real axis except for one simple zero on real axis 
at M = 0. The Bethe Ansatz equations follow from (14.81) . 

a{vk + —) = -l, k = l,... ,M = N + u-l. (4.9) 
Following the steps of section 3, we remove the root of T[u) at the origin by defining 

f(«) = ^, (4.10) 

where //(m) is a function whose only real root is a simple zero at the origin, with /i(0) = 
, /i'(0) 7^ 0. In terms of T(m), the T — Q equation becomes 

T{u) = Uu) ^ ^(1 + aiu)) = t+(n) ^(1 + aiu)) , (4.11) 



where 



^ sinh(2n±z/.)sinh(n±f)sinh(.Tf)^(,) ^ ij^ 

^{u) 2 
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Exploiting the analyticity of \nT{u) near the real axis, we have 



du [lnT(u)]"e 



II iku 



C 



C2 



l" Aku 



ill iku 




/ du [ln(l + a(u))]"e 



" Aku 



e*'^"+ / du [ln(l + a(u))]"e 

C2 



" Aku 



(4.13) 



where the contour C is again chosen as in Figure 1. 

We first carefully separate the extra-roots terms in Q{u) before proceeding with 

further computation involving Fourier transforms defined in ( ]3.13p . Introducing ( ]4.6p in 
f lO]) . we get 



1 

Q{'^) = W T sinh(u — Wk) sinh(u + Wk) 



k=l 

v-l 



X 



X 



k=l 



smh [-{u -wl' - —)) smh [-{u + wl' + —)) 



(1) , ^vr. 



sinh {^{u - w^}^ + ^)) sinh {^{u + 



w 



(1) 



)) 



(4.14) 



Defining 



2 ^ 



(4.15) 



k=l 



and using fl3.13p and f l3.26p . we arrive at the following result for the Fourier transform of 
[lng(«)]", 



v-1 
2 



LQ"{k) = LQ"{k) + -^^Y^cosH—) cos{kwf^] 



(4.16) 



sinhfvrfc) ' 2 

j=i 

Using the periodicity for Q{u), namely 

Q{u) = Q{u — in), u E Ci, and Q{u + ijj) = Q{u + iji — in), m G C2 (4.17) 
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and the identity fl3.26p . the second and the fifth terms in fl4.13p reduce to 



du < In 



Ci 



Q{u + ifi{iy - 1)) 



e*'^" + i du {\n 

'C2 



Q{u — ifilu — 1)) 



iku 



-Ae-^ cosh(^) sinh (-(vr - fi))LQ"{k) 



8nk 



sinh(7r/c) 



smh{kfi) cosh(— ) cos(A;u']"'^-*^ 



As in section 3, we define 



C{k)= / du[\nU{u)]"e""'+ / du[\nt4u)]" e 



" Aku 



(4.18) 



(4.19) 



As a result, we obtain the following from fl4.13p . 



LQ"{k) 



7rk 

e 2 



+ 



4cosh(^f)sinh((7r-^)|) 
47rA;sinh(^)e^ 



LA"{k) + CA"{k) + C{k) 



kiT , 



cosh( — ) cos(A; 



sinh(7rA;) sinh ((vr — /i)|) ' 2 

Following definitions have again been adopted, 

A{u) = 1 + a{u) , A{u) = 1 + a{u) 



w 



(4.20) 



(4.21) 



Next, we proceed to derive the NLIE for the lattice sine-Gordon model. Fourier transform 
of a{u) (see f l4.7p ) along C2 yields. 



La"{k) 











/„**■ 









Q{u + iiJ,{v - 1)) 
_Q{u - z/i(z/ - 1)) 



Aku 



sinh(2M + in) sinh(M + f ) sinh(M - f^) B^+\u) 0(m + f ) 



sinh(2M - i/i) sinh(M - f ) sinh(M + ^) ^{-)(m) . 



Aku 



(4.22) 



17 



Using fHTT]) and IK2^ . we obtain 
L^'{k) = LQ"{k)[e-^''' -e^''-^^''] 

v-l 



— — — — sinh(/cu) cosh( — ) cos 
^ ^ 7=1 



{k w 



+ du {\n 

'C2 



sinh(2M + ifi) sinh(u + f ) sinh(M - fi) + f ) 

sinh(2M - ifi) sinh(M - f ) sinh(M + fi) 0(m - f ) 



Aku 



(4.23) 



Inserting fl4.20p into (14.231) yields the following NLIE for the sine-Gordon model in Fourier 
space, 

L^'{k) = G{k) \LA"{k) + CA"{k) \ + Cr^k) 



Airk 



v-l 
2 



cosh(^) sinh (|(7r - /i)) 



sinh(^) cosh(— ) cos(A; 



where 
G(A;) 

Cr(fc) 



sinh ((tt - 2/i)|) 
2cosh(f )sinh ((7r-/i)|) ' 

G(A:)C(A:) 



(4.24) 



(4.25) 



+ du < In 



C2 



'sinh(2u + z^) sinh(M + f ) sinh(M - ^) i3(+)(M) 0(m + 
sinh(2M - ifi) sinh(u - f ) sinh(M + fi) B(~){u) 



Aku 



(4.26) 



Adopting similar steps of the previous section, and using the following for C{k) (which 
is obtained from fl4.19p ). 



Cik) 



du < In 



sinh(2M + tfi) sinh(M + f ) sinh(M - ^) B^+\u) 0(m + f ) 

u}0{u-'-§) 



C2 1^ L^inh(2M - ifi) sinh(M - f) smh{u + ^) i3(-)( 
du [ln/i(M)]"e'^" 



Aku 



(4.27) 
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we find, 

. _,^JiVcos(M)_smh(f)cosh((2,-.)|) 
^ ' \ cosh(f) cosh(f )sinh((^-7r)|) 

cosh(^) sinh ((2/i - 7r)|) sinh cos{kiJ,b) 



cosh(f ) sinh ((/i - 7r)|) cosh(f ) sinli ((/i - 7r)|) 

+ geosh(^) cos(fc^.f )| . (4.28) 

cosh(f )sinh((7r-/i)|) ^ ^2^ ^ ^ ^ 

The second term in f l4.27p reduces to —271/2. Converting ( I4.24p to coordinate space and 
integrating twice, we obtain 

/CO roo 
du' G{u — u' + ie) ln(l + a{u — it)) — I du' G{u — u — ie) ln(l + d{u' + ie)) 
■00 J —00 



sinh 



TVU 



where we recall that G{u) is the Fourier transform of G{k) (see f l4.25p ) 

1 ^ 

G{u) = — dk e-*'^" Gik) , (4.30) 

and Pbdryiu) is given by 

Pbdryiu) = du'R{u') = - du'Riu'), (4.31) 

where /?(«) refers to the Fourier transform of R{k) given below, 

sinh(^) cosh ((2/i -7r)|) 



R{k) = -27r 



cosh(^) sinh ((/i — 7r)|) 
cosh(^) sinh ((2/i - 7r)|) sinh (|/i) cos{kfib) 



cosh(f ) sinh ((/i - 7r)|) cosh(f ) sinh {{fi - 7r)|) 



2 sinh(^) ^ fcvr 

+ m Vf^ -xvV cosh(— )cos(A;wfM. (4.32) 

cosh(f )sinh((7r-/i)|) ^ ^2^ ^ ^ ^ 

The factor in f l4.29p is obtained as before by considering the m — )■ oo limit of ( 14. 7p and 
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Next, taking the continuum limit (A— j-oOjA^— J-oOjA— i-O) and defining the renormal- 
ized rapidity 6 as in f l3.34p . f l4.29p becomes 



Ina(^) 



dO' G{e -e' + le) ln(l + a{e' - ie)) - / dO' G{e -e' - is) ln(l + a(^' + is)) 



i2mL sinh + i PhdryiG) + 



ITT 



U-1 



(4.33) 



where definitions f l3.36p have been used. By making following identification 

fie) ^ H-m) 



(4.34) 



fl4.33p takes the following form, 

/oo 
de' '^mG{e-e' - ie) ln(l - e^(^'+*")) 
'OO 

ZTT 

- iPbdryiO) + iTT - — -y . 



2imL sinh 6 



(4.35) 



4.2 Boundary and Casimir energies 

Now we compute the boundary correction (order 1) and the Casimir correction (order 1/L) 
to the energy. Using (13. 390 - 03. 410 and (13.430 . where now 



T^^^(m) = sinh(2'u ± i/i) sinh(M ± — ) sinh('u =1= — —) B''^\u) (f){u ± — 



X 



Q{u ± i/i(z/ — 1)) 

we have the following for the energy E, 

g d 



(4.36) 



E 



, , lnT(+)(n + ^)-lnT(-)(n-^^ 
Adu \ ^ 2' ^ 2' 



g f dk 



-ikA 



A / 2n 



'e'^LTM'ik)-e-'^LT(^'ik) 



(4.37) 



where now one finds, 



fj.k 



e^LT(+y{k) - e-2 LT(-)'(/t) 



j^Sw'e^) + (e^'^-"^' - e-^'')I4'{k) + (e^'^-i)^ - e"^'^)^^^ 

{-ik) 



e-'^li^'ik) + Ae-^ sinh((7r - fi)^)LQ'{k) 



2 



IQm , T /^^N /; (1)n 



smh[Tik) 



sinh( — ) cosh( — ) cos(A; 



(4.38) 
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LQ'{k) is determined from fl4.20p and LQ'{k) = jpr^LQ"{k). As a result, after some algebra 
fl4.38p becomes, 



e'^LT(+y{k) - e-^LT(-)'(A;) 



2tanli(^)e~^ cosh ((/x - ^)k)L(p{k) 



+ tanh(^) 



2'Ki ^ 27rz cosh((yU — j)k) 



sinh ^ 



sinli ^ 



+ 



+ 



47ri cosli((/i — 



sinh ^ 



+ e-^LB(-y{k) + e^LBM'ik) 



coshf 



LA'(A;) + £A'(A;) 



(4.39) 



Note that the terms with extra-roots {wj^''} cancel out eventually. Introducing (I4.39p in 
(I4.37p . one has the following result for the energy, 

E = El + Ei + EyL (4.40) 

where 

El = 



r dk e-*'=^tanh(^) he""? cosh((/i - -)k)L^{k) 
271 A J_^ 2 L 2 . 



(4.41) 



El 



9 



27rA 



dk e 



-ikA 



tanh(^) 



27ii 27ri cosh((/i — |)fc) 



sinh 



jik 



+ 



2 

Ani cosh((/i — j)k) 
sinh ^ 



sinh ^ 



+ e 



■3H'(A;) + e'^3(+)'(A;) 

(4.42) 



2nA 



dk e 



-ikA 



coshf 



LA'{k) + £A'(A;) 



(4.43) 



Equations (I4.4ip . fl4.42p and fl4.43p represent the bulk, boundary and Casimir energies re- 
spectively. The bulk energy E^, is computed as before (refer to section 3.2) to yield 



1. 2 
Er = -Lm cot — 

4 2/i 



(4.44) 



where L and m are defined by (12. 4p . We now consider the boundary energy (14.420 . Using 
results fl33T]) (with B^^^) and ([S^SD, we have 



e'"^ m-y (k) + e"^ m+y (k) 



— ; — - cosh( — ) cos(/e/xo) 
sinh Ilk 2 



(4.45) 
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Substituting (14.451) into (14.421) . we obtain the boundary energy 

(4.46) 



2 



vrz/ 7TU cosh (bn) 
1 + cot h2cot 2- ^ ' 



4 ■ 2 sin^ 



where we have used the same contour as for the bulk energy, namely by closing the contour in 
the lower half plane and selecting only the contribution from the pole at A; = — ^. Comparing 
(I4.46P with Al. Zamolodchikov's result [TT| [T2] for the energy of the continuum sine-Gordon 
model (single boundary), 



777 



I cos (7r/(2A)) + I sin (7r/(2A)) - ^ + cos(r//A) + cosh(i9/A) 



2cos(7r/(2A)) [2 ^ ' ^ 2 ^ ' ^ 2 

and using A = ^ — 1 = we conclude that 

= A7r6 = (A + l)/3 , (4.48) 

One can verify that (I4.48P coincides with the result in [T9] with a± = 1 (refer to equation 
(3.16) of that paper) for the constraint case (see equations (12.81) and (I2.10p ) studied in section 
3 of this paper. This suggests that the relation between the boundary parameters of the 
continuum model {ri± , and the lattice model {a± , b±) in [T9] might hold true in general. 
Consequently, the relation between the boundary parameters of the lattice model {a± , f3±) 
and that of the continuum action (/i-t ,(p^) (see (12. 3p ) as given in (12. 5p should also be true 
for the case with arbitrary I3± , hence possibly indicating that ( 12. 5p might hold true in general 
as well. One notes that for i/ = 3 , 5 , 7 , . . ., the boundary energy becomes 

m 



(4.47) 



El 



1 + (-1) — (1 - 2cosh(67r)) . (4.49) 



Lastly, we consider the Casimir energy (14. 43 p . Passing to coordinate space and taking 
the continuum limit, we obtain 

Ec = ^ I du^m\ — — — ^ — 1 ln(l + a(u + ^e)) (4.50) 

cosh ^(A — M — xe) 




which reduces to 



Ec = / c/M5>meM("+"hn(l + a(M + 2e)) (4.51) 

2/i J -00 

in A — 7- 00 limit (where (^ cosh^(A-M-ie) ) ^ ^e^'i^^^^^^'"'^^) and after using (12. 4p . Similar 
manipulations of section 3.2 reduces Ec to (We have used the fact that A(m) = A(— m) and 
'^mz = —Qmz in the process.) 

Ec = -— d9Qmsmh{9 + ie)\n{l + a{9 + ie)) (4.52) 
27r 
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Invoking the identification f{6) = ln(— 3(6*)), the above becomes 

Ec = / de smh{9 + is) ln(l - e-^(^+^")) . (4.53) 

By considering the UV hmit where L — )• and utihzing the steps in [9], we obtain the 
following for the Casimir energy Ec, 



Ec=-^^ (4,54) 



where the central charge c is given by 



1-^ P^r, 00 +vr 



= l-6(^) (4.55) 

where Pbdryi^o) = ~^(2 + ^;rj-)' following section, we present numerical results for the 

central charge of the corresponding open spin-| XXZ spin chain, which agree with f l4.55p . 
We emphasize that for the case considered in section 3, the central charge of the sine-Gordon 
model (in the UV limit, L — )■ 0) was found to coincide with that of the XXZ spin chain. See 
section 3.2.1 in [19] for more details. 



4.3 Numerical results 

In this section, we present some numerical results and the extrapolated value for the central 
charge of XXZ spin chain. Such results are obtained as follows: We first numerically solve 
the Bethe equations (14.91) for some finite number of spins. We use the solution to calculate 
Casimir energy from the following 

E Ei)ulk ~\~ Eijoy^Tfifidyy -\- E(j(igiYnir (4.56) 

In (14561) . E is given by [24j 

1 1 1 

E = -sinVV ■■ ^ + -(iV- l)cosu (4.57) 

2 ^^cosh(t;,-f)cosh(t;fc + f) 2^ ^ ^ ^ ^ 



Having determined the Bethe roots numerically, one uses known expressions for Ebuik [38 
namely 



r-oo 1 1 

2 ' - J- -L 



Ebulk = -iVsin> / dX- — — — —— + -Ncosfx, (4.58) 

, [cosh(2yuA) — cosyuj cosh(7rA) 2 
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and Ei 



boundary 



boundary 



given by 

sin /i 



du 



+ 



J~oo 2cosh(a;/2) 

sinh(u;/2) cosh((z/ - 2)a;/2) 
sinh(i^a;/2) 



+ 



sinh((i^ - 2)uj/4) 
2sinh(z/a;/4) 

smh.{u/2) cos{bu) 



-cos/i. (4.59) 
2 



to determine Ecasimir- Then using -Ecasimir 
central charge, Cg// for that value of A^, 



IT Sin fi 



sinh(z/a;/2) 

one determines the effective 



24:iJ.N ^'^ff-' 



-eff 



{E-E, 



bulk 



Et 



boundary ) 



(4.60) 

TT^ sm fi ' 

Since we are ultimately interested in — )■ oo limit, we employ an algorithm due to Vanden 
Broeck and Schwartz [33]- [31] to extrapolate these values for central charge to — )■ oo 
limit. Table 1 below shows the Cg// values for few finite even A^, computed for some values 
of fi and b, for the lowest energy state that we considered in section 4.1 of this paper. The 
extrapolated values (-8.02719 and -7.97641) obtained from the Vanden Broeck and Schwartz 
algorithm agree with the result obtained (= -8) from f l4.55p . 



A^ 


Ceff, b = 1.65 


Ceff, b = 1.79 


112 


-3.0586758038329473 


-2.510905382065773 


120 


-3.1514592759743896 


-2.6001056254878296 


128 


-3.2385273078233214 


-2.6842360822480313 


136 


-3.3204955969641814 


-2.7638181748627413 


144 


-3.3978856780639934 


-2.8392958325047744 


152 


-3.4711433285531212 


-2.911050348946888 


160 


-3.5406526531820894 


-2.9794118013959547 


168 


-3.6067470178385475 


-3.0446679508432566 


176 


-3.6697176469097252 


-3.1070712683288697 


184 


-3.7298204634077208 


-3.1668445456417547 


192 


-3.787281586552167 


-3.2241854218624653 


extrapolated value 


-8.02719 


-7.97641 



Table 1: Central charge Cg//, for z/ = 3 for two different boundary parameter values from 
numerical computations based on A^ = 112 ,120 ,. . . ,192 and extrapolated values at A^ — )■ oo 
limit (Vanden Broeck and Schwartz algorithm). 

The agreement between the calculated and the extrapolated values indicates that as with 
the constraint case studied in [H], the central charge of the sine-Gordon model (in UV limit). 
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coincides with that of the corresponding spin-| XXZ quantum spin chain. Furthermore, 
the numerical results also indicate that the extrapolated value of Cg// is independent of 
the boundary parameter, as expected for models with Neumann boundary condition. This 
implies that the results for Cg// and the conformal dimension A = ^ ^^^^ = ^^J^-^^ have more 
resemblance to spin chains with diagonal boundary terms, as one would expect from the 
dependance [3 HH 113 SS] rather than [20] which is the anticipated form for the 
conformal dimensions for spin chains with nondiagonal boundary terms. For more complete 
and detailed discussion on this, readers are urged to refer to |l5l |20] . 

5 Discussion 

From the proposed T — Q equation of an open XXZ quantum spin chain with nondiagonal 
boundary terms, we derived the NLIE for the lowest energy state of an open spin-i XXZ/sine- 
Gordon model with two boundaries. We first rederived the NLIE for the case where the lattice 
boundary parameters obey certain constraint, which was treated in [19l |20]. However, in 
contrast to the approach used there, we employed a method that utilizes the T — Q equation 
of the spin chain model. We next derived the NLIE for the case without such a constraint 
among the boundary parameters, where one of the lattice parameters is set to be completely 
arbitrary. The lowest energy state of this spin chain model has complex sea of Bethe roots 
which is rather a feature common to critical spin-1 XXZ spin chain. We also obtained 
the boundary energy and Casimir energy for the lowest energy state. We then presented 
relations between the boundary parameters of the continuum model and that of the lattice 
model which coincide with the ones found in [19] for the constraint case, hence suggesting 
that these relations might hold true in general. 

Having found the 1/L correction, we proceeded to compute the central charge of the 
sine- Gordon model in the UV limit. We also solved the corresponding spin-| XXZ chain 
numerically for some finite values of A^. We used the solutions to compute correction 
for these values, then extrapolate them to the N ^ oo limit using Vanden Broeck and 
Schwartz algorithm. The extrapolated value of the effective central charge, Cg// is found to 
be in agreement with that of the sine-Gordon model in UV limit. 

In addition, the numerical results also indicate that Cg// is independent of the boundary 
parameters, as expected for models with Neumann boundary condition. The result for the 
conformal dimension A, turns out to be similar to that of the XXZ spin chain models with 
diagonal boundary terms rather than the nondiagonal ones, to which the model studied here 
belongs to. Such a feature however had been encountered before in literature [12] • There it 
was pointed out that such a behaviour can be possibly attributed to spectral equivalences 
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between diagonal-nondiagonal open XXZ spin chains [39]-|41]. 

There are other problems that one can explore and address further. For example, one 
could investigate the open spin-1 XXZ chain with general integrable boundary conditions in 
a similar manner, along the line of |31J, since solutions for such a model are already available 
[13| m] . It will also be interesting to study boundary excitations for these cases. We look 
forward to address these issues in near future. 

Acknowledgments 

I am grateful to the referee for pointing out a crucial reference. 

References 

[1] M. Gaudin, "Boundary Energy of a Bose Gas in One Dimension," Phys. Rev. A4, 386 
(1971); La fonction d'onde de Bethe (Masson, 1983). 

[2] F.C. Alcaraz, M.N. Barber, M.T. Batchelor, R.J. Baxter and G.R.W. Quispel, "Surface 
exponents of the quantum XXZ, Ashkin- Teller and Potts models," J. Phys. A20, 6397 
(1987). 

[3] C.J. Hamer, G.R.W. Quispel and M.T. Batchelor, "Conformal Anomaly And Surface 
Energy For Potts And Ashkin-Teller Quantum Chains," J. Phys. A20, 5677 (1987). 

[4] E.K. Sklyanin, "Boundary conditions for integrable quantum systems," J. Phys. A21, 
2375 (1988). 

[5] H.J. de Vega and A. Gonzalez-Ruiz, "Boundary K-matrices for the six vertex and the 
n{2n - 1) An-i vertex models," J. Phys. A26, L519 (1993) | hep-th/9211114 . 

[6] S. Ghoshal and A.B. Zamolodchikov, "Boundary S-Matrix and Boundary State in Two- 
Dimensional Integrable Quantum Field Theory," Int. J. Mod. Phys. A9, 3841 (1994) 
[hep-th/9306002] . 

[7] P. Fendley and H. Saleur, "Deriving boundary S matrices," Nud. Phys. B428, 681 
(1994) |hep-th/9402045 . 

[8] M. T. Grisaru, L. Mezincescu and R. I. Nepomechie, "Direct calculation of the boundary 
S matrix for the open Heisenberg chain," J. Phys. A28, 1027 (1995) |hep-th/9407089j 



26 



[9] A. LeClair, G. Mussardo, H. Saleur and S. Skorik, "Boundary energy and bound- 
ary states in integrable quantum field theories," Nucl. Phys. B453, 581 (1995) 
|hep-th/9503227 . 

[10] A. Doikou and R. I. Nepomechie, "Direct calculation of breather S matrices," J. Phys. 
A32, 3663 (1999) | hep-th/9903066] 

[11] Al. Zamolodchikov, invited talk at the 4th Bologna Workshop, June 1999. 

[12] Z. Bajnok, L. Palla and G. Takacs, "Finite size effects in boundary sine-Gordon theory," 
Nucl. Phys. B622, 565 (2002) |hep-th/0108157] . 

[13] J.-S. Caux, H. Saleur and F. Siano, "The Josephson current in Luttinger liquid- 
superconductor junctions," Phys. Rev. Lett. 88 (2002) 106402 fcond-mat/0109103l . 

[14] T. Lee and C. Rim, "Thermodynamic Bethe Ansatz for boundary sine-Gordon model," 
Nucl. Phys. B672, 487 (2003) |hep-th/0301075] . 

[15] J. -S. Caux, H. Saleur and F. Siano, "The two-boundary sine-Gordon model," Nucl. 
Phys. B672, 411 (2003) c ond-mat/6 306328|. 

[16] C. Ahn, M. Bellacosa and F. Ravanini, "Excited states NLIE for sine-Gordon 
model in a strip with Dirichlet boundary conditions," Phys. Lett. B595, 537 (2004) 
|hep-th/03 12176 . 

[17] J. Cao, H.-Q. Lin, K.-J. Shi and Y. Wang, "Exact solutions and elementary excitations 
in the XXZ spin chain with unparallel boundary fields," [cond-mat/0212163]; 
J. Cao, H.-Q. Lin, K.-J. Shi and Y. Wang, "Exact solution of XXZ spin chain with 
unparallel boundary fields," Nucl. Phys. B663, 487 (2003). 

[18] R.L Nepomechie, "Functional relations and Bethe Ansatz for the XXZ chain," J. Stat. 
Phys. Ill, 1363 (2003) |hep-th/02110 01]; 

R.L Nepomechie, "Bethe Ansatz solution of the open XXZ chain with nondiagonal 
boundary terms," J. Phys. AST, 433 (2004) |hep-th/0304092|; 

R.L Nepomechie and F. Ravanini, "Completeness of the Bethe Ansatz solution of the 
open XXZ chain with nondiagonal boundary terms," J. Phys. A36, 11391 (2003); Ad- 
dendum, J. Phys. AST, 1945 (2004) |hep-th/0307095] . 

[19] C. Ahn and R.L Nepomechie, "Finite size effects in the XXZ and sine-Gordon models 
with two boundaries," Nucl. Phys. B6T6, 637 (2004) [hep-th/0309261j . 



27 



[20] C. Ahn, Z. Bajnok, R.I. Nepomechie, L. Palla and G. Takacs, "NLIE for hole excited 
states in the sine-Gordon model with two boundaries," Nucl. Phys. B714, 307 (2005) 
|hep-th/0501047 . 



[21] A. Kliimper, M.T. Batchelor and P.A. Pearce, "Central charges of the 6- and 19-vertex 
models with twisted boundary conditions," J. Phys. A24, 3111 (1991). 

[22] C. Destri and H. de Vega, "New thermodynamic Bethe Ansatz equations without 
strings," Phys. Rev. Lett. 69, 2313 (1992) |hep-th/9203064] ; 

C. Destri and H. de Vega, "Unified approach to thermodjTiamic Bethe Ansatz and fi- 
nite size corrections for lattice models and field theories," Nucl. Phys. B438, 413 (1995) 
|hep-th/9407117l; 

C. Destri and H. de Vega, "Non-linear integral equation and excited-states scaling func- 
tions in the sine-Gordon model," Nucl. Phys. B504, 621 (1997) ^hep-th/9701107] . 

[23] D. Fioravanti, A. Mariottini, E. Quattrini and F. Ravanini, "Excited state Destri-De 
Vega equation for sine-Gordon and restricted sine-Gordon models," Phys. Lett. B390, 
243 (1997) |hep-th/960809l]; 

G. Feverati, F. Ravanini and G. Takacs, "Nonlinear integral equation and finite volume 
spectrum of sine-Gordon theory," Nucl. Phys. B540, 543 (1999) | hep-th /9805117); 
G. Feverati, "Finite volume spectrum of sine-Gordon model and its 



restrictions," hep-th/0001172 



[24] R. Murgan and R.I. Nepomechie, "Bethe ansatz derived from the functional rela- 
tions of the open XXZ chain for new special cases," J. Stat. Mech. P05007, (2005) 
[hep-th/0504124|; Addendum, J. Stat. Mech. P11004, (2005) |hep-th/0504124]; 
R. Murgan, R.I. Nepomechie and C. Shi, "Boundary energy of the open XXZ chain 
from new exact solutions," Ann. Henri PoincareT, 1429 (2006) [hep-th/0512058j . 

[25] P. Di Vecchia and S. Ferrara, "Classical solutions in two-dimensional supersymmetric 
field theories," Nucl. Phys. B130, 93 (1977); 

J. Hruby, "On the supersymmetric sine-Gordon model and a two-dimensional 'bag' 
Nucl. Phys. B131, 275 (1977). 

[26] S. Ferrara, L. Girardello and S. Sciuto, "An infinite set of conservation laws of the su- 
persymmetric sine-Gordon theory," Phys. Lett. B76, 303 (1978); 

L. Girardello and S. Sciuto, "Inverse scattering like problem for supersymmetric mod- 
els," Phys. Lett. B77, 267 (1978); 

R. Sasaki and I. Yamanaka, "Supervirasoro algebra and solvable supersymmetric quan- 
tum field theories," Prog. Theor. Phys. 79, 1167 (1988). 



28 



[27] R. Shankar and E. Witten, "The S matrix of the super symmetric nonhnear sigma 
model," Phys. Rev. D17, 2134 (1978). 

[28] C. Ahn, D. Bernard and A. LeClair, "Fractional supersymmetries in perturbed coset 
CFTs and integrable soliton theory," Nucl. Phys. B346, 409 (1990); 
C. Ahn, "Complete S matrices of supersymmetric sine-Gordon theory and perturbed 
superconformal minimal model," Nucl. Phys. B354, 57 (1991). 

[29] C. Dunning, "Finite size effects and the supersymmetric sine-Gordon models," J. Phys. 
A36, 5463 (2003) hep -th/02 10225] ; 

Z. Bajnok, C. Dunning, L. Palla, G. Takacs and F. Wagner, "SUSY sine-Gordon theory 
as a perturbed conformal field theory and finite size effects," Nucl. Phys. B679, 521 
(2004) |hep -th/0369120 . 

[30] A. Hegedus, F. Ravanini and J. Suzuki, "Exact finite size spectrum in super sine-Gordon 
model," Nucl. Phys. B763, 330 (2007) |hep-t h/0610 012]. 

[31] C. Ahn, R.I. Nepomechie and J. Suzuki, "Finite size effects in the spin-1 XXZ and 
supersymmetric sine-Gordon models with Dirichlet boundary conditions," Nucl. Phys. 
B767, 250 (2007) |hep -th/0611136] . 

[32] J. Suzuki, "Spinous in magnetic chains of arbitrary spins at finite temperatures," J. 
Phys. A32, 2341 (1999); 

J. Suzuki, "Excited states nonlinear integral equations for an integrable anisotropic spin 
1 chain," J. Phys. A37, 11957 (2004) |hep-th/0410243]. 

[33] J-M van den Broeck and L.W Schwartz, "One-parameter family of sequence transfor- 
mations," SIAM J. Math. Anal. 10, 658 (1979) 

[34] C.J. Hamer and M.N. Barber, "Finite-lattice extrapolations for and Z5 models," J. 
Phys. A14, 2009 (1981) 

[35] M. Henkel and G. Schiitz, "Finite-lattice extrapolation algorithms," J. Phys. A21, 2617 
(1988) 

[36] Al.B. Zamolodchikov, "Mass scale in the sine-Gordon model and its reductions," Int. 
J. Mod. Phys. AlO, 1125 (1995). 

[37] N.Yu. Reshetikhin and H. Saleur, "Lattice regularization of massive and massless inte- 
grable field theories," Nucl. Phys. B419, 507 (1994) [he p-th/9309135j . 



29 



[38] C.N. Yang and CP. Yang, "One-Dimensional Chain of Anisotropic Spin-Spin Interac- 
tions II. Properties of the Ground-State Energy Per Lattice Site for an Infinite System," 
Phys. Rev. 150, 327 (1966). 

[39] J. de Gier, A. Nichols, P. Pyatov and V. Rittenberg, "Magic in the spectra of the XXZ 
quantum chain with boundaries at A = and A = —1/2," Nucl. Phys. B729, 387 
(2005) [ hep-th/0505062j . 

[40] A. Nichols, V. Rittenberg and J. de Gier, "One-boundary Temperley-Lieb algebras in 
the XXZ and loop models," J. Stat. Mech. P03003, (2005) licond-mat/ 041 1512 . 

[41] Z. Bajnok, "Equivalences between spin models induced by defects," J. Stat. Mech. 
P06010, (2006) [hep-th/0601107j . 

[42] R. Murgan, "Finite-size correction and bulk hole-excitations for special case of an open 
XXZ chain with nondiagonal boundary terms at roots of unity," JHEP 05, 069 (2007) 
[hep-th/0704.2265]. 

[43] L. Frappat, R.I. Nepomechie and E. Ragoucy, "Complete Bethe ansatz solution of 
the open spin-s XXZ chain with general integrable boundary terms," J. Stat. Mech. 
P09008, (2007) [math-ph/0707.0653]. 

[44] R. Murgan, "Bethe ansatz of the open spin-s XXZ chain with nondiagonal boundary 
terms," JHEP 04., 076 (2009) [hep-th/0901 . 3558]. 

[45] H. Saleur, "Lectures on non pertubative field theory and quantum impurity problems," 
[cond-mat/9812110j. 

[46] I. Affleck, M. Oshikawa and H. Saleur, "Quantum brownian motion on a triangular 
lattice and c = 2 boundary conformal field theory," Nucl. Phys. B594, 535 (2001) 
lcond-mat/000 90841 . 



30 



